Strain-induced interface reconstruction in epitaxial heterostructures 
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We investigate in the framework of Landau theory the distortion of the strain fields at the interface 
of two dissimilar ferroelastic oxides that undergo a structural cubic-to-tetragonal phase transition. 
Simple analytical solutions are derived for the dilatational and the deviatoric strains that are globally 
valid over the whole of the heterostructure. The solutions reveal that the dilatational strain exhibits 
compression close to the interface which may in turn affect the electronic properties in that region. 

PACS numbers: 71.27.+a, 81.30.Kf 



I. INTRODUCTION 

Recent discoveries in material science related to sev- 
eral unexpected properties of epitaxial heterostructures 
made of different transition metal oxide (TMO) materi- 
als, bring in the forefront of interest the problem of inter- 
face reconstruction through the developement of sponta- 
neous strain at the interfacei*^. Lattice distortion close 
to the interface is known to result in charge redistribu- 
tion that leads to the formation of a two-dimensional 
electron gas (2DEG) and metallicity in that region^^. 
Most of the TMOs of interest are ferroelastics that un- 
dergo structural transitions^ from a cubic/pseudocubic 
to a lower symmetry phase with decreasing tempera- 
ture. Notably, heterostructures containing strontium ti- 
tanate (SrTiOa), a band-insulator oxide undergoing a 
cubic-to-tetragonal (CTT) structural transition at Tg ~ 
105 ^K, exhibit extraordinary interfacial properties be- 
low Ts] metallicit}*ii^i^i^, superconductivity^, and non- 
linear Hall effect Moreover, in LaTiOs/SrTiOa and 
LaAlOs/SrTiOa heterostructures, the structural transi- 
tion of SrTiOa causes the overlayers to stabilize in a 
tetragonal phase with an in-plane lattice constant almost 
equal to that of SrTiOa close to the interfac e -^-^i-^^ . 

It has been discussed in the past that the electro- 
magnetic properties of TMOs couple to the elastic de- 
grees of freedon>i^— . The effect of tensile and com- 
pressive strains to the electronic conduction properties 
at the interface of TMO heterostructures has already 
been addressed experiment all}*i^iii. Furthermore, strong 
polarization enhancement in ferroelectric TMO super- 
lattices driven by interfacial strain has been unambi- 
giously observed^^. In the present work we apply con- 
tinuous elasticity theory through a Ginzburg-Landau de- 
scription in terms of the strain tensor components to 
a heterostructure. Based solely on symmetry consid- 
erations, Ginzburg-Landau theory can provide a reli- 
able description of the equilibrium behavior of a system 
near a phase transition. It has been recently used to 
show theoretically the emergence of a multiferroic state 
of a EuTiOa film on (LaA103)o.29(SrAli/2Tai/203)o.7i 
(LSAT) substrate^^, to provide a physical understanding 



of the strain-induced metal-insulator phase coexistence in 
manganites^^, and to explain phase separation between 
metallic ferromagnetic and insulating charge-modulated 
phases^S. 

We investigate the interfacial effects on the strain- 
state of a bilayer heterostructure, composed of dissimilar 
TMOs that join at a single planar interface and propose 
a strain-based mechanism that may help understand the 
formation of a 2DEG. In particular, we obtain approx- 
imate analytical solutions for the dilatational and the 
deviatoric strain fields in the bilayer, that exhibit spatial 
variation due to breaking of the uniformity. Notably, the 
dilatational strain field exhibits a well-defined minimum 
at the interface corresponding to local compression^. We 
argue that the suppression of the dilatational strain field 
in the interfacial region may encourage the formation of 
a 2DEG. The proposed strain-based mechanism does not 
exclude other possible mechanisms, like, e.g., the orbital 
and/or the electronic reconstruction mechanisms— i^^. 



II. GINZBURG-LANDAU THEORY AND 
EQUATIONS OF MOTION 

In the Lagrangian description of elasticity the sym- 
metric strain tensor is defined as e^j = {^i.j + /2 
= x^y^z), where Uij is the j—th. derivative of the 
i—th component of the displacement vector u of a ma- 
terial point relative to its position in the parent phase. 
The six symmetry adapted strains for the CTT structural 
transition are defined as^^ 

63 = ^ ^y^y ~ ^^^'^)' ^4 = ^{Uy,z + Uz^y), (2) 

while 65 and cq are given by 64 with cyclic permutation 
of the indices. The deviatoric strains 62 and 63 form 
the two-component order parameter (OP) of the CTT 
transition. Both the OP and the non-OP strains are 
coordinate-independent in the uniform product (tetrago- 
nal) phase in static equilibrium, with the latter custom- 



arily being set to zero. In a TMO heterostructure, where 
the uniformity of the product phase is broken due to the 
interface, ah e^'s vary spatiahy; in that case, their second 
derivatives are linked through compatibihty relations^-. 
In a non-uniform state, the non-OP strains cannot be all 
set to zero. Specifically, in TMO heterostructures the di- 
latational strain ei, which is concommitant to 63^^^^^, ex- 
hibits measurable compression indicating its importance 
in their structural properties^. 

In ferroelasticity theory, the strain energy density T 
of a material undergoing a CTT structural transition is 
expanded in powers of the invariants of the strain ten- 
sor and their products around the energy of the parent 
phase^^^^. Thus, the functional T is expressed solely in 
terms of the e^'s and their spatial derivatives. Guided by 
previous works we adopt a functional T of the form 
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where the Ginzburg-Landau coefficients 
^1,^2,^4,^7, and 61,62,63 are related to the second-, 
third-, and fourth-order elastic coefficients of the parent 
phase, respectively, through (in Voigt notation) 
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= 3(Cii+2Ci2) 
= 2(Cii - C12), 



while di , d2 are two independent strain-gradient coeffi- 
cients. In accordance with common principles of Lan- 
dau theory, the critical temperature dependence of the 
62 elastic constant, 62 ex (T — Tg), is supposed to be true 
close to the transition point. 

In a single material at static equilibrium, the spatially 
homogeneous strains in the product phase are the lowest 
energy solutions of the conditions dJ^/dci = 0. Neglect- 
ing the non-OP terms, the energy density landscape on 
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FIG. 1: (color online) Strain energy density lanscape on 
the 62 — 63 plane for = 25.4, 64 = 225, and C2 = —10 
(left); —50 (right), exhibiting the familiar pattern of three 
degenerate minima. 



the 62 — 63 plane exhibits the familiar pattern of three 
degenerate minima corresponding to three different vari- 
ants in the tetragonal phase (Fig. 1). Notably, a non- 
zero 61 preserves the energetic degeneracy of the three 
variants. We are particularly interested in the variant 
having 62 = 0; this is because TMO heterostructures are 
usually grown along the z— direction and both materials 
go into a 6— tetragonal phase at low temperatures^^^^^. 
For 62 = 0, the energy landscape on the 61 — 63 plane 
shown in Fig. 2 exhibits significant qualitative differences 
for different a2 values. Specifically, the strain 61 varies 
from positive (i.e., expansional) to negative (i.e., com- 
pressional) with decreasing the magnitude of a2. A small 
a2 absolute value is however expected from the princi- 
ple 61 << 62,63 for all martensitic transformations'^!^. 
Moreover, a small a2 leads to negative 61 and 63, in ac- 
cordance with the empirical principle for the ferroelastic 
transitions of close-packed solids, i.e., that cooling of the 
solid is usually accompanied by a decrease of volume. 
The dependence of the strains on 62 for all three variants 
is shown in Fig. 3 for two different values of a2. We 
later refer to the two materials forming the bilayer het- 
erostructure, which occupy the regions z < and z > 0, 
as the left (L) and the right (R) material, respectively. 
The Ginzburg-Landau parameters used in Figs. 1-3 are 
those given for the left material in Table I, and they 
have been calculated from the corresponding elastic co- 
efficients through Eqs. (jH). The second- and third-order 
elastic coefficients of the left material are those reported 
for SrTi03^^^^^, while for the fourth-order ones a reason- 
able choise was made (Table I). Note that for the param- 
eter a2, which can be treated as a phenomenological one, 
we have also used values that are smaller than the one 
given in Table I for the left material (i.e., a2 = —31). 

The dynamics of the displacements is governed by the 
Euler-Lagrange equations 
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FIG. 2: (color online) Strain energy density lanscape on the 
ei — 63 plane for C2 = —50, and a2 — —31 (left); —10 (right). 
The other Ginzburg-Landau parameters are those for the left 
material given in Table 1. 



where 



on 



(6) 



are the strain tensor and the dissipative strain tensor, 
respectively, po is the density in the parent phase, and 1Z 
is the the Rayleigh dissipation function 



^ = ^^1^1 + ^^2(^2 + el) + y,{el + + ej). (7) 
Then, from Eq. ^ we get 

poili = di^i + ^{csHi + c'^Hi), (8) 

where 
and 

Wi ciei + — Ue2 + ' eici - —63(11) 

with q = +1 (—1) for i = x {y). The functions W/ 
and Gi have the same form with that of the W^'s, with 
the obvious change Cj and Cj dj, respectivelly 

{j = 1,2), and the i^^'s are lengthy nonlinear functions 
of ei, 62, and 63, which are given in Appendix A. 



III. APPROXIMATE SOLUTIONS AND 
INTERFACE RECONSTRUCTION 

In order to separate the interfacial effects on the strain- 
state of the bilayer heterostructure (from those origi- 
nating from external boundaries, domain walls, disloca- 
tions, etc.), we consider two monodomain, semi-infinite 
TMOs joined along a chemically abrupt, planar interface 
at z = 0. Eqs. (|8|) could be simplified in a strict way. 
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FIG. 3: (color online) The strains ei (black-solid), 62 (red- 
dashed), 63 (green-dotted), as a function of C2 for the three 
variants, for a2 = —7 (left panels); —31 (right panels). The 
other Ginzburg-Landau parameters are those for the left ma- 
terial given in Table I. 



since at low temperatures the strains depend on one co- 
ordinate only and 62 = 0^^'^^. However, for non-zero 
ei the simpification of Eqs. (jS]) following the strict way 
is a non-trivial task, which makes preferable the use of 
a simple ansatz for the displacements. This ansatz as- 
sures that the strains depend only on the z— coordinate 
and that 62, as well as the small strains 64, 65, ee, are 
identically zero. Assume that the strains exhibit a rela- 
tively strong z— coordinate dependence in the proximity 
of the interface, while they attain their static equilibium 
values for large enough \z\. This approximation seems 
well-suited for heterostuctures composed of TMOs with 
small lattice mismatch (i.e., LaTi03/SrTi03). Indeed, 
both experimental observations^ (discussed below) and 
first-principles calculations^ indicate that strain inhomo- 
geneity and lattice deformation occur within a few layers 
near the interface. Thus, for practical purposes, it is suf- 
ficient for the two layers of the heterostructure to be thick 
enough for the deformation to vanish relatively far from 
the interface. The choise of semi-infinite layers was made 
only for mathematical convinience. 
We then introduce the ansatz 



2^' 



hz^f{z), (12) 



where f{z) is a yet unknown function, and 

<^ = -^(eio + V^eso), 6 = +-^(eio - 2^8630), (13) 

with eio and 630 being the values of ei and 63, respec- 
tively, far from the interface. The non-zero strains are 
then 

ei = eio + /(^), 63 = 630 - -^/(^), (14) 

where the prime denotes differentiation with respect to 
z. Substitution of Eq. (p!l|) into Eqs. (|8|) results, in the 
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static limit, in the equation 



di + |)G"'=(ci + |)G' + i?;(G), (15) 



where G = G{z) = f'{z)^ and Rz with 62 = is 
o 2 , T 3 , / ^2 \ 2 , / ^7 ^4 \ 3 
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where ei and 63 are meant to be expressed in terms of 
/ through Eq. ([14]). After rearrangement, Eq. (fT5]) 
becomes 



G'" = G'(/^ + 3AG + 6/iG^), 
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FIG. 4: (color online) The strains ei and 63 as a function of ; 



for C2 



-l,c,^ 



-0.76 (black-solid), C2 



-10, 



-7.56 



(red-dashed); C2 = —50, = —37.8 (green-dotted), and (a)- 
(b): = -31, = -29.5; (c)-(d): = -10, = -9.5. 
The other Ginzburg-Landau parameters are those for the left 
and the right material given in Table I. 



where 



3di 
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with 



ctiz = ci^^C2 = Cii, diz = di^^d2, (19) 
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Eq. (p!7|) can be reduced to a quadrature that has the 
analytic solution 

Q± = 4^ ^g±V^(^T^o) ^ ^gTv^(^T^o) ^ ^ (21) 

where —A = 4/^/i — A^, Ai = —2 A, and zq is a constant 
of integration. Integration of gives 

f± = tanh-1 , ^ ^ + C±, (22) 

where G^ are constants of integration. 

The displacements and the strains in each material of 
the bilayer can be written in terms of and G^ from 
Eqs. (fT2]) and (p!4|) . respectively. Specifically, the solu- 
tions for ei, 63 in the material left, (right) from the 
interface occupying the region z < {z > 0) are written 
as 



MR) 



.T^G^^-\ e3^(-^=e3^o^-^--^G+(-)(24) 



where the superscript L {R) indicates the value of the 
corresponding quantity in the left (right) material. For 
this choice, the integration constants in Eq. (|22]) are 

G^ = tanh~^ (^ 2^R.fi ) ' ^^^^ ^ derivatives 
vanish on either side of the heterostructure far from the 
interface, in accordance with our earlier assumptions. 

In order to obtain solutions for Uz, ei, that are glob- 
ally valid over the whole bilayer structure, we impose the 
following (internal) boundary conditions at the interface 



(25) 



where the stress component azz is obtained from ^z^ Eq. 
([To]), in the static limit, and z'^ is the location of the 
interface that is not necessarily at zero. We thus distin- 
guish between the positions of the actual interface, where 
the strains exhibit significant variation, and the interface 
which is the natural boundary of the two materials. The 
actual and the natural interfaces could be slightly dis- 
placed one another due to reconstruction of the interface, 
similarly to that observed in Ag(lll)/Ru(0001)^^. Eqs. 
([25|) can be satisfied for appropriate values of zq and z^ 
which can be obtained numerically. 

The strains ei and 63 along the z— direction, that is 
perpendicular to the interface, are shown in Fig. 4 for 
several combinations of C2 and a2. The strain ei exhibits 
a minimum close to the interface, indicating relative lat- 
tice compression in that region. Notably, compressionally 
strained layers at the PbTi03/SrTi03 interface, corre- 
sponding to a reduced c— axis lattice parameter of the 
PbTi03 film in the first few unit cells, have been ex- 
perimentally observed^^. This interface reconstruction is 
solely due to the elastic properties of the materials of the 
bilayer. The dependence of ei on 2: and C2 is shown in 
the left panel of Fig. 5. The corresponding dependence 
of the well's depth, and the constants z^ and zq, is 
shown in the right panel of Fig. 5. Thus, with decreasing 
C2 (i.e., becoming more negative) D increases, while zq^ 
which is a measure of the well's width, decreases. Also, 
decreases with decreasing C2, so that the actual in- 
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FIG. 5: (color online) Left panel: The strain ei as a function 
of the 2;— coordinate and the C2 parameter. Right panel: The 
depth of the well D, and the numerically obtained constants 
zq and as a function of C2. The other Ginzburg-Landau 
parameters are those for the left and the right material given 
in Table I. 



terface approaches the natural one at low temperatures. 

The values of the Ginzburg-Landau parameters used 
in the calculations of the strains in Figs. 4 and 5 are 
given in Table I, calculated from the corresponding elas- 
tic coefficients through Eqs. For the left (L) mate- 
rial, the second- and third-order constants are those re- 
ported for SrTiO s^-^1^^ , while the second-order constants 
for the right {R) material are those reported recently for 
LaTiOs^^. The values of the elastic coefficients, whose 
values are not reported in the literature due to the lack of 
experimental and theoretical data, are chosen to be rea- 
sonable for perovskites (Table I). However, the obtained 
interfacial effects persists for a wide range of values of 
the higher-order elastic coefficients. The strain gradient 
coefficients di and d2 are treated as phenomenological pa- 
rameters to set the scale of the deformed region. Their 
values for the left and the right material are chosen to 
be, respectively, = ^2 = 6 and = d2 = 7, in 
units of 10 ~^ N. For this choise, significant variations 
in ei and 63 occur within ~ 1 nm corresponding to '^2-3 
TMO layers. The value of C2 was taken to be negative 
for both materials, in accordance with common practices 
of Ginzburg-Landau theory of phase transitions. In Figs. 
4 and 5, where C2 and/or a2 vary for both materials, a 
constant ratio between the left and the right value was 
assumed, calculated from their values given in Table I. 
The parameter a2 is again treated phenomenologically, so 
that pairs of a2 and a2 with absolute values smaller that 
those given in Table I (but with the same ratio) have been 
used. These values are also consistent with the principle 
ei << 62,63 for all martensitic tranformations. 

It has been reported that LaA103 and LaTi03 follow 
the structure of the SrTi03 substrate when the latter un- 
dergoes a GTT transition^^^^^. The correlation of these 
structural changes and the electromagnetic properties in 
these systems may be empirically seen in the observation 
of an enhancement in the interfacial charge carrier mo- 
bility and magnetization below Tj^^-'^. In TMO het- 
erostructures undergoing a GTT transition, significant 



lattice deformation occurs at the interface region due 
to lattice mismatch, that results in spontaneous strains. 
Despite the empirical evidence for the effects of interfa- 
cial lattice deformation on the electromagnetic properties 
of TMO heterostructures, the relationship between the 
strain and the formation of a 2DEG remains largely un- 
explored. In a Ginzurg-Landau approach that includes 
charge and/or magnetic degrees of freedom, the dilata- 
tional strain 61 couples linearly to the charge densit}^. 
Then, the results of Fig. 4 reveal that 61 serves as an 
effective potential well which may affect the charge dis- 
tribution throughout the heterostructure. In particular, 
the compressed interfacial region may attract and confine 
electron charges. The localized charges may contribute 
to the formation of a 2DEG, a prerequisite for interfacial 
metallicity in TMO heterostructures. 
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225 


393 
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250.0 


244.0 




b7 


-35.4 


-40.3 



TABLE I: Second-, third-, and fourth-order elastic coefficients 
(in units of 10^^ N/m'^), and Ginzburg-Landau (GL) dimen- 
sionless coefficients for the left and the right material used in 
the calculations. 



IV. CONCLUSIONS 

We applied continuum elasticity to investigate theoret- 
ically the strain-state of bilayer TMO heterostructures 
within a Landau theory, and we have obtained simple 
approximate solutions for the fields 61 and 63. Interface 
reconstruction may lead to electronic charge redistribu- 
tion in the heterostructure, and particularly to electronic 
charge concentration in the interface region favoring the 
formation of a 2DEG. The presence of a minimum in 
the dilatational strain field demonstrates that possibil- 
ity, linking thus the elastic to the electronic properties of 
TMOs. Although such a reconstruction is a microscopic 
phenomenon involving significant changes of atomic ar- 
rangements at the interface^^, it results in macroscopic 
changes of the unit cells that can be observed experimen- 
tally. Those changes can be described, at least qualita- 
tively, by the Ginzburg-Landau theory, and their impli- 
cations on the electron charge distribution of the bilayer 
can be inferred from basic physical laws. 
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Appendix A: Nonlinear functions 



The nonlinear functions Ri {i = x^y, z) are given by 
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It can be easily checked that for 62 = we have Rx = Ry 
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